Calculated spin-orbit splitting of all diamond-like and zinc-blende semiconductors: 
Effects of pi/2 local orbitals and chemical trends 



^1- 
O 

o 

< 
O 

CO 



CZ3 



c 

o 

CN 
> 
On 
O 

m 
O 

o 

-I— > 
c3 



i 

c 

o 
o 



X 



Pierre Carrier 1 and Su-Huai Wei 1 

1 National Renewable Energy Laboratory, Golden CO 80401, USA 
(Dated: February 2, 2008) 

We have calculated the spin-orbit (SO) splitting Aso = e(r"8u) — e(Tr v ) for all diamond-like group 
IV and zinc-blende group III-V, II- VI, and I- VII semiconductors using the full potential linearized 
augmented plane wave method within the local density approximation. The SO coupling is included 
using the second variation procedure, including the pi/ 2 local orbitals. The calculated SO splittings 
are in very good agreement with available experimental data. The corrections due to the inclusion of 
the pi/2 local orbital are negligible for lighter atoms, but can be as large as ~250 meV for 6p anions. 
We find that (i) the SO splittings increase monotonically when anion atomic number increases; (ii) 
the SO splittings increase with the cation atomic number when the compound is more covalent such 
as in most III-V compounds; (iii) the SO splittings decrease with the cation atomic number when 
the compound is more ionic, such as in II- VI and the Ill-nitride compounds; (iv) the common-anion 
rule, which states that the variation of Aso is small for common-anion systems, is usually obeyed, 
especially for ionic systems, but can break down if the compounds contain second-row elements such 
as BSb; (v) for IB- VII compounds, the Aso is small and in many cases negative and it does not 
follow the rules discussed above. These trends are explained in terms of atomic SO splitting, volume 
deformation-induced charge renormalization, and cation-anion p-d couplings. 

PACS numbers: 71.55.-i, 61.72. Vv; 78.20.Bh, 78.40.-q 



I. INTRODUCTION 

Spin-orbit (SO) splitting Aso = ^8v) — £{^rv) at 
the top of the valence band of a semiconductor is an 
important parameter for the determination of optical 
transitions in these systemsi 2 * 3 . It is also an impor- 
tant parameter to gauge the chemical environment and 
bonding of a semiconductor pLiiLi Extensive studies 
of SO splitting, both theoretically 8 ! 9 ! 10 ! 11 ! 12 ! 13 : 14 ! 15 and 
experimentally, 16 . 17 ! 18 ! 19 ! 20 . 21 ! 22 ! 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 ! 30 have 
been carried out in the past. However, most of these 
studies focussed on a specific compound or a small group 
of similar compounds. Therefore, the general trends of 
the spin-orbit splitting in zinc-blende semiconductors is 
not very well established. From the experimental point 
of view, some of the data were measured more than 30 
years ago^i 7 - and the accuracy of these data is still under 
debate. For example, previous experimental data suggest 
that CdTe and HgTe have SO splittings Aso at about 
0.8 and 1.08 eV, respectively^ These values have been 
used widely by experimental groups^ to interpret opti- 
cal and magneto-optical transition data of CdTe, HgTe, 
and related alloys and heterostructures. However, re- 
cent experimental data suggest that Aso for CdTe and 
HgTe are instead around 0.95 eVS 7 . and 0.91 eV^& With- 
out basic understanding of the general trends of varia- 
tion of Aso m tetrahedral semiconductors, it is difficult 
to judge what should be the correct value of Aso f° r 
CdTe and HgTe. There are also several non-conventional 
II- VI and III-V semiconductors that do not have a zinc- 
blende ground state (e.g., CdO, MgO, GaBi, InBi), but 
that do form zinc-blende alloys with other compounds, 
and are currently under intensive research as novel opto- 
electronic materials ^M^iSl Therefore, it is important 



to know the spin-orbit splittings of these compounds in 
the zinc-blende phase and understand how they vary as 
a function of alloy concentration x in the alloy. 

From the theoretical point of view, various approxi- 
mations have been used to calculate and/or predict SO 
splitting Aso- However, it is not clear how these ap- 
proximations affect the calculated Aso- For example, 
one of the most widely used procedures for calculating 
the SO coupling using the density functional theorjiS 5 . 
(DFT) and local density approximation^ 7 . (LDA) is the 
second-variation method 38,3 — used in many all-electron 
linearized augmented plane wave (LAPW) codesi 4 ^* 41 * 4 ^ 
In this approach, following the suggestion of Koelling and 
Harmon, 38 the Hamiltonian of the relativistic Dirac equa- 
tion is separated into a "J-weighted-averaged" scalar rel- 
ativistic Hamiltonian Hsr, in which the dependancy on 
the quantum number k [where k — ±(j + 1/2), with 
\j\ = \l + s\ = Z =p 1/2] is removed from the full Hamilto- 
nian, and a spin-orbit Hamiltonian Hso with 
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where 



M = m - 



e-V 



is the relativistically enhanced electron mass, c is the 
speed of light, V is the effective potential, e is the eigen- 
value, and s and I are the Pauli spin and angular mo- 
mentum operators, respectively. The scalar relativis- 
tic Hamiltonian, which includes the mass velocity and 
Darwin corrections, is solved first using standard diag- 
onalization method for each spin orientation (or solved 
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FIG. 1: Comparison of P1/2, P3/2, an d Pi=i orbitals in atomic 
As and Bi showing the large discrepancy between p\ji and the 
p !=1 orbitals, especially for the heavier Bi atom. 



just once if the system is not spin polarized). The SO 
Hamiltonian is included subsequently, in where the full 
Hamiltonian is solved using the scalar relativistic wave- 
functions as basis set. Normally, only a small number 
of scalar relativistic wavefunctions are included in the 
second step, and only the spherical part of the poten- 
tial within a muffin-tin sphere centered on each atomic 
site is used in the SO Hamiltonian. The advantage of 
the second- variation method is the physical transparency 
(e.g., it keeps spin as a good quantum number as long as 
possible) and the efficiency, because, in most cases in the 
second step, only a small number of basis functions are 
needed to have good agreement with solutions of fully rel- 
ativistic Dirac equations. This approach has been shown 
to obtain Aso that is in excellent agreement with ex- 
periments. For example, the calculated Aso for GaAs is 
0.34 eV compared with experimental data of 0.34 eVU 
However, one major approximation in the "J-weighted- 
averaged" treatment is the replacement of the two P1/2 
and P3/2 orbitals by one pi = i orbital. Although this is a 
good approximation for atoms with low atomic number, 
it has been show that such approximation fails for heavy 
atoms i^iMiiii The main reason for this failure is because 
the P1/2 orbital has finite magnitude at the nuclear site, 
whereas the I = 1 orbital has zero magnitude at the nu- 
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FIG. 2: Chemical trend of the spin-orbit splittings for all 
diamond-like group IV and zinc-blende group III-V, II- VI, 
and I- VII semiconductors, including the p 1( / 2 local orbitals. 
The graph corresponds to the data in column "LAPW+P1/2" 

of Tables hi irn and nm 



clear site. Figure plots the P1/2, P3/2, and pi = \ orbitals 
for As (Z=33) and Bi (Z=83). As we can see, thep!/ 2 or- 
bital deviates significantly from the pi = i orbital near the 
origin. The error clearly increases as the atomic number 
increases, and is very large for heavier elements such as 
Bi. Therefore, the P1/2 orbital is not very well repre- 
sented near the nuclear site using the pi = \ orbital, even 
with the addition of its energy derivative in the lineariza- 
tion procedure^ Consequently, the SO splitting cannot 
be accurately evaluated, in general, with solely the pi = \ 
orbital. However, no systematic studies have been done 
to evaluate the effect of the pi / 2 orbital on the calculated 
SO splitting A50 
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FIG. 3: Atomic spin-orbit splittings £(^3/2) — ^(Pi/2) f° r 
atoms studied in this paper. The spin-orbit splittings increase 
as a function of the atomic number Z. See Table ITVl for data 
subdivided according to their respective groups. 

The objective of this paper is to do a systematic study 
of the SO splitting Aso of all diamond-group IV and 
zinc-blende groups III-V, II- VI, and I- VII semiconductors 
using the first-principles band-structure method within 
the density functional formalism. We find that the cal- 
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culated SO splittings including the pi/2 local orbital are 
in good agreement with available experimental data. The 
general chemical trends of the Aso are revealed and ex- 
plained in terms of atomic SO splittings, volume effects, 
and p-d coupling effects. 



II. METHOD OF CALCULATIONS 

The calculations are performed using the full potential 
linearized augmented plane wave (FLAPW) method as 
implemented in the WIEN2k code42*^ The frozen core 
projector augmented wave (PAW) approach implemented 
in the VASP cod o 4 ^ 4 ^ is used for comparison. We used 
the Monkhost-Pack 4 ^ 4 x 4 x 4 k points for the Bril- 
louin zone integration. For the FLAPW method, SO 
coupling is included using the second-variation method 
performed with or without the pi/2 local orbitals. Highly 
converged cutoff parameters in terms of the numbers of 
spherical harmonics inside the muffin-tin region and the 
plane waves in the interstitial region, as well as local 
orbitals for low-lying valence band states (anion s and 
cation d states), are used to ensure the full convergence 
of the calculated values. For the PAW method, high pre- 
cision energy cutoffs have been chosen for all semiconduc- 
tors (as large as 37 Rydberg for the nitrides and oxides). 

In most cases, the band structure calculations are per- 
formed at the experimental lattice constants. For com- 
pounds that have only experimental lattice constant in 
the wurtzite structure, such as ZnO, we assume that 
zinc-blende ZnO has the same volume as in its wurtzite 
structured For BSb, the (Al, Ga, In)Bi, and the (Be, 
Mg, Cd, Hg)0, that does not have either zinc-blende or 
wurtzite experimental structure parameters, the LDA- 
calculated lattice constants are used. For the silver 
halides and the gold halides, the LDA lattice constants 
have been corrected according to the small discrepancy 
between LDA and experiment values of Agl (more pre- 
cisely, 0.088 A have been added to the LDA lattice con- 
stants of silver halides and gold halides). The LDA- 
calculated lattice constants are expected to be reliable. 
For example, our predicted^ lattice constant of GaBi is 
a = 6.324 A, whereas recent experimental observation^ 
finds a value around 6.33±0.06 A, in good agreement 
with our prediction. All the lattice constants used in our 
calculation are listed in Tables HI ITU and lllll 



III. EFFECT OF THE p 1/2 LOCAL ORBITAL 

Tables UJ [n] and IIIII present the calculated SO split- 
tings data for all diamond-like group IV and zinc-blende 
groups III-V, II- VI, and I- VII semiconductors. The cal- 
culated values are obtained with or without the p\/2 local 
orbitals. We find that including the p 1 / 2 local orbital pro- 
vides a better variation basis for the T^ v state, lowers the 
eigen energy, and, therefore, increases the SO splitting 



TABLE I: Calculated spin-orbit splitting Aso for all diamond 
group IV and zinc-blende group III-V semiconductors, using 
the FLAPW method with or without the P1/2 local orbitals 
and the frozen-core PAW method. Our results are compared 
with available experimental data. Our error analysis suggests 
that the uncertainty of the LDA calculated value is less than 
20 meV. 
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Aso — e (r8u) - e (^7v)- The correction due to the pi/2 
orbital increases as the atomic number increases. Since 
the VBM consists of mostly anion p state, the dependence 
is more on anion atomic numbers. We find that correc- 
tions due to the inclusion of the P1/2 local orbital (for 
both anions and cations) are negligible for lighter atoms, 
are ^10 meV for 4p anions, ~40 meV for bp anions and 
can be as large as ~250 meV for 6p anions. Thus, for Bi 
compounds (AIBi, GaBi, and InBi), large errors could be 
introduced if the pi/2 local orbital is not included^ In 
all the cases, inclusion of pi/2 local orbital brings a bet- 
ter agreement between the calculated Aso and available 
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FIG. 4: Charge distribution at the VBM for SiC. The charges 
are mostly distributed on the carbon atom site. 

experimental data. 



IV. CHEMICAL TRENDS 

Figure [3 shows the general chemical trends of the cal- 
culated SO splittings Aso for all diamond-like group IV 
and zinc-blende III-V, II- VI, and I- VII semiconductors, 
with inclusion of the P\j2 local orbitals. We find that 
(i) the SO splittings increase monotonically when anion 
atomic number increases; (ii) the SO splittings increase 
with the cation atomic number when the compound is 
more covalent, such as in most III-V compounds; (iii) 
the SO splittings decrease with the cation atomic num- 
ber when the compound is more ionic, such as in II- VI 
and the Ill-nitride compounds; (iv) for compounds with 
the same principal quantum number, A50 increases as 
the ionicity of the compounds increases. Finally, (v) 
the halides (IB- VII) constitute a special case because the 
VBM in IB- VII is no longer an anion p dominant stated 
Therefore, IB- VII compounds do not follow the rules dis- 
cussed above. 

To understand these chemical trends, we will first dis- 
cuss the factors that can affect the SO splitting A50 
for the systems studied here, (a) Dependence on the 
atomic number. The atomic SO splitting between the 
P3/2 and P1/2 states increases as a function of atomic 
number Z. Table Hvl gives the calculated splitting of the 
atomic fine structures, £(^3/2) - e(Pi/2)i as a function of 
the atomic number Z in their respective groups. Fig- 
ure |31 (related to Table IIV() shows the variation of the 
atomic spin-orbit splittings as a function of the atomic 
numbers, for all atoms considered. The spin-orbit split- 
tings increase with the atomic number, as expected. 4 ^ 
The increases follow approximately a power law with 



TABLE II: Calculated spin-orbit splitting A so for all IIA-VI 
and IIB-VI semiconductors, using the FLAPW method, with 
or without the p\/2 local orbitals, and the frozen-core PAW 
method. The lattice constants with a * corresponds to one at 
their LDA energy minimum (for ZnO**, the lattice constant 
of the zinc-blende structure is chosen so that its volume is 
equal to that in the wurtzite structure). Our results are com- 
pared with available experimental data. Our error analysis 
suggests that due to the overestimation of the p-d hybridiza- 
tion, our calculated Aso is underestimated by 30, 40, and 
110 meV for Zn, Cd, and Hg compounds, respectively. For 
other compounds, the LDA error is estimated to be less than 
20 meV. 
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5.0162* 
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364 


369 


370 
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6.4820 
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810 c , 800 d , 900 f 


HgO 


5.1566* 


-285 


-281 


-292 




HgS 


5.8500 


-100 


-87 


-108 




HgSe 


6.0850 


235 


254 


238 


450 c , 396 d 300 9 


HgTe 


6.4603 


762 


800 


781 


1080 c , 910 s 



a Reference 22! 



b Reference 
c Reference I 
d Reference 
e Reference 
f Reference 
9 Reference 



Aso(P3/2 — P1/2) Z a , where a is close to 2. (b) De- 
pendence on the volume. As the volume of the com- 
pound decreases, the charge distribution in the crystal is 
renormalized. The bonds become more covalent. More 
charge is pushed into a region near the nuclei. Because 
the SO coupling is larger near the nuclear site, the SO 
splitting Aso usually increases as the volume decreases, 
(c) Dependence on the cation valence d orbital. 
The VBM in a majority of zinc-blende semiconductors 
consists of mostly anion p and a smaller amount of cation 
p orbitals. By symmetry, the VBM state in zinc-blende 
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(b) BSb 




FIG. 5: Charge density of the VBM state for GaSb and BSb, 
showing that for BSb the role of cation and anion is reversed. 

structure can couple with the cation tid orbitals. The 
cation t^d orbital has a negative contribution^^ to the 
SO splitting Aso (i-e., the T Sv is below the Ty v state). 
Thus, large mixing of heavy cation d orbitals in VBM 
can reduce Aso- 

Using the discussion above, we can now understand 
the general chemical trends of the SO splitting Aso- 

(i) The SO splittings increase monotonically when an- 
ion atomic number increases. For example, Aso in- 
creases from 13 — > 49 — > 302 — > 697 meV when the atomic 
number increases from C — > Si — > Ge — > a-Sn; from 12 
— > 86 — + 342 — > 738 -^2150 meV when the anion atomic 
number increases from GaN — > GaP — » GaAs — > GaSb — > 
GaBi; from -60 -> 50 -> 369 -> 880 meV when the an- 
ion atomic number increases from CdO — ► CdS — ► CdSe 
— > CdTe; from -85 — > 82 — > 455 when the anion atomic 
number increases from CuCl — > CuBr — > Cul. This is 



TABLE III: Calculated spin-orbit splitting Aso for all IB- VII 
compounds, using the FLAPW method, with or without the 
Pi/2 local orbitals, and the frozen-core PAW method. Our 
results are compared with available experimental data. We 
use experimental lattice constants— for CuX (X= CI, 
Br, I) and Agl. The lattice constants for the other AgX and 
AuX compounds are estimated from calculated LDA lattice 
constans and experimental lattice constant of Agl. Due to 
the overestimation of the d character in the VBM, the LDA 
underestimate the Aso by 20, 60, and 170 meV for chlorides, 
bromides, and iodides, respectively. 
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a References 54' andl55t 



because the VBM has large anion p character, and the 
atomic SO splitting of the anion valence p state increases 
with the atomic number (see Table UVj l. One of the inter- 
esting case is SiC. The calculated Aso of 14 meV for SiC 
is very close to the one of diamond (13 meV), indicating 
that SiC is a very ionic material with its VBM containing 
mostly C character. Figure 21 depicts the contour plot of 
the charge distribution at the VBM for SiC, which shows 
that the VBM charge is located on the carbon atom site. 

(ii) The SO splittings increase with the cation atomic 
number when the compound is more covalent, such as in 
most III-V compounds. For example, Aso increases from 
216 — > 300 — > 342 — > 352 meV when the atomic number 
increases from BAs — > AlAs— > GaAs — > InAs; from 366 — > 
681 — > 738 — > 755 meV when the atomic number increases 
from BSb -> AlSb -> GaSb -» InSb. This is because 
for covalent III-V compounds, the VBM contains signif- 
icant amount of cation p orbitals. Therefore, when the 
cation atomic number increases, the SO splitting Aso 
also increases. It is interesting to note that A so f° r BX 
(X=P, As, and Sb) is significantly smaller than that for 
their corresponding common-anion compounds. For ex- 
ample, Aso(BSb)=366 meV is only about half of the 
value of Aso(GaSb)=738 meV. This is because boron 
is much more electronegative than other group III ele- 
ments. Thus, BX compounds are much more covalent 
than the other III-V semiconductors. Figure compares 
the charge distribution of the VBM states for BSb and 
GaSb. We see that for GaSb, most of the VBM charge is 
on Sb atom site, whereas for BSb, a large portion of the 
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FIG. 6: Charge density of the VBM states for Ge, GaAs, 
ZnSe, and CuBr showing that as ionicity increases, the charge 
is more localized on the anion site. For ZnSe and CuBr, it 
also shows antibonding d character on the Zn and Cu sites, 
respectively^ 

VBM charge is on the B atom site. Because boron has 
a small atomic number (Z=5), the SO splitting of B 2p 
states is very small, leading to very small Ago for BX. 
This indicates that the common-anion rule, which states 
that the variation of Ago is small for common anion 
systems, does not apply to all BX, which are extremely 
covalent. 

(iii) The SO splittings decrease with the cation atomic 
number when the compound is more ionic, such as in 
II- VI and Ill-nitride compounds. For example, Ago de- 
creases from 449 — > 399 meV when the atomic number 
increases from BcSc — » MgSe; from 965 — > 869 meV when 
the atomic number increases from BeTe — + MgTe; from 
398 — > 369 — > 254 meV when the atomic number in- 
creases from ZnSe — > CdSe — > HgSe; from 21 — > 19 — > 
12^0 meV when the atomic number increases from BN 
to A1N -> GaN -> InN. This is because for ionic II- VI 
and Ill-nitride systems, the VBM is mostly an anion p 
state, thus the Ago is not sensitive to the cation atomic 
number or potential. However, when cation atomic num- 
ber decreases, say from Mg to Be, the volume of the 
compounds decreases (Table |nj, and therefore, due to 
the charge renormalization effect, the Ago increases. In 
particular, for the IIB-VI systems and the Ill-nitrides, 
the coupling between cation d and anion p also plays an 
important role in the observed trend, because the p-d hy- 
bridization is significant in these systems (See Fig. [BJ;) . 
The p-d hybridization reduces Asor^ and the effect 
increases when cation atomic number increases. This ex- 
plains why Ago(HgX) (for X=0, S, Se, Te) is smaller 
than Ago(CdX), even though they have similar volume, 



TABLE IV: Atomic SO splitting e(p 3 / 2 ) - ^{Pi/2) for the com- 
pounds of Tables [I] EI an d Hill according to their atomic 
groups. The data are also depicted in Fig. [31 as a function of 
atomic numbers Z. 



element 


atomic number 
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Sn 


Z=50 


463 


V 
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282 


Sb 


Z=51 


632 


Bi 


Z=83 


1 968 
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Z= 8 


37 


s 


Z=16 


86 


Se 


Z=34 


386 


Te 


Z=52 


815 


VII 






CI 


Z=17 


127 


Br 


Z=35 


509 


I 


Z=53 


1 029 



and why Ago (InN) is smaller than Ago (GaN). Note that 
negative Aso can exist in some of the compounds such 
as ZnO, CdO, and HgO where the anion is light, so their 
p orbitals have only a small contribution to Ago, but the 
negative contribution of the cation d orbital is large. 

(iv) For compounds with the same principal quantum 
number n, Aso increases as the ionicity of the compound 
increases. For example, for n = 2, from C — > BN — > 
BeO, the SO splittings Ago increase from 13 — > 21 — > 
36 meV; for n = 3, from Si -> A1P -> MgS, the SO 
splittings increase from 49 — * 59 — * 87 meV; for n = 4, 



from Ge — ► GaAs — » ZnSe, the SO splittings increase 
from 302 -> 342 -> to 398 meV; for n = 5, from a-Sn 
— > InSb — > CdTe, the SO splittings increase from 697 
— > 755 — > 880 meV. The reason for this increase can be 



understood from plots in Fig. which show the charge 
distribution of the VBM states of Ge, GaAs, and ZnSe. 
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As the system changes from group IV — > III-V — > II- VI, 
the compound becomes more ionic and the VBM becomes 
more localized on the anion site with increasing atomic 
number, thus Ago increases. It is interesting to note 
that the differences of Ago between the II- VI, the III-V, 
and the group IV compounds in the same row increases 
as n increases (almost doubles when n increases by one) . 
This is explained by the fact that the atomic number Z 
almost doubles when n is increased by one, whereas the 
atomic SO splitting is proportional to Z a with a close 
to 2 (See Table llVl and the discussion above); thus, the 
difference is proportional to Z. 

(v) The A IB X vn halides (A IB = Cu, Ag, Au; X vn = 
CI, Br, I) constitute a group of special compounds that 
do not follow the rules discussed above. For example, 
when moving from ZnSe to CuBr with increased ionicity 
(see Figure EJ, the SO splitting of CuBr (82 meV) is 
much smaller than that for ZnSe (398 meV). The SO 
splitting of Agl (664 meV) is also much smaller than that 
of CdTe (880 meV). Furthermore, many of the IB- VII 
compounds (CuCl, AgCl, AuCl, and CuBr) have negative 
SO splittings, and for these ionic compounds CuX vn has 
much smaller SO splittings than AgX vn and AuX vn . 
The origin of these anomalies is due to the fact that for 
most of the IB- VII compounds the VBM is no longer 
an anion p dominated state. Instead, they are cation d 
states strongly hybridized with the anion p state. For 
instance, in Figure Efd) we show that the VBM of CuBr 
has a very pronounced antibonding d character at the 
cation Cu site. Because the d state has negative Ago, 
this explains why some of the IB- VII compounds have 
negative Ago- Furthermore, because Cu 3d level is much 
higher than Ag 4d and Au 5d levels, the VBM of Cu 
halides contains more cation d character than Ag and 
Au compounds. This explains why Cu halides have much 
smaller Ago than the Ag and Au common anion halides. 



V. COMPARISON WITH EXPERIMENTS 



understand the origin of the discrepancy, we performed 
the following tests. First, we considered a different 
numerical approach, i.e., the frozen core PAW method 
as implemented in the VASP code to calculate the SO 
splitting Ago. Despite the large difference in the way 
the SO coupling is implemented in the calculations, we 
find that the Ago calculated with the PAW method 
are very similar to that obtained with the FLAPW 
method. For a-Sn and HgTe, the results obtained by 
the PAW method are 689 and 781 meV, respectively, in 
good agreement with the FLAPW-calculated values of 
697 and 800 meV. Next, we estimated the effect of p-d 
coupling. It has been argued that the LDA-calculated 
cation d orbitals are too shallow^ so p-d hybridization 
at the VBM is overestimated, which may lead to smaller 
calculated Ago- To verify if this is the possible reason, 
we performed the following calculations. (i) After 
obtaining the converged LDA potential, we removed the 
cation d orbital from the basis set to calculate the Ago- 
We find that for a-Sn, this procedure has no effect on 
the calculated Ago- This is consistent with the fact that 
for this compound, the cation d and anion p separation 
is large enough that the amount of cation d orbital 
at the VBM is not sufficient to affect the calculated 
Ago- For ZnTe, CdTe, and HgTe, removing the cation 
d orbital increases the Ago by 48, 63, and 253 meV, 
respectively. These values are the upper limit on the 
possible effect of p-d coupling on the calculated Ago- 
(ii) To get more reliable estimates on the LDA error of 
the calculated Ago, we added an external potential^ 
on the cation muffin-tin sphere to push down the cation 
d orbitals such that the calculated cation binding energy 
is close to the experimental photemission data^ In this 
case, the calculated Ago is 0.94, 0.91, and 0.90 eV for 
ZnTe, CdTe, and HgTe, respectively. The above analysis 
demonstrates that the possible LDA error in calculating 
Ago is less than 30, 40, and 110 meV for Zn, Cd, and 
Hg compounds, respectively, and much smaller for other 
compounds. 



Our calculated results with the p 1 / 2 lo- 
cal orbitals are compared with experimental 

^.[■^ 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27.28. 29. 30 p Qr most 

semiconductors the agreement is very good. For ex- 
ample, the calculated value for diamond (13 meV) is 
in very good agreement with the recent experimentally 
derived value of 13 meV^S The experimental value 
for SiC in the zinc-blende structure (10 meV)iiii is 
smaller than the one for C, therefore, does not follow the 
chemical trend. We suggest that the measured value is 
possibly underestimated. For most semiconductors, the 
difference between theory and experiment is usually less 
than 20 meV. However, there are several noticeable cases 
in which the difference is much larger. For example, 
for a-Sn, the calculated value is 697 meV, whereas the 
value in experiment data™ is ~800 meV. For HgTe 
the calculated value at 800 meV is much smaller than 
the widely used experimental valued of 1080 meV. To 



Our analysis above suggests that Ago for a-Sn and 
HgTe should be around 0.70 and 0.90 eV, respectively, 
smaller than the experimental values of 0.80 and 1.08 
eV, respectively. The origin of this discrepancy is still 
not very clear. But we notice that a-Sn and HgTe 
are semimetals, i.e., the T^ c state is below the VBM. 
This makes the accurate measurement of the Ago for 
these compounds more challenging. Indeed, recent 
measurements^ of Ago for HgTe show that it has a value 
of 0.9 eV, in good agreement with our predicted value. 
We also notice that the recent reported experimental SO 
splitting for InSb)2i which has a very small band gap 
(0.24 eV), agrees well with our calculation. Further ex- 
perimental studies are needed to clarify these issues. 
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VI. SUMMARY 

In summary, we have studied systematically the SO 
splitting A so of all diamond-like group IV and zinc- 
blende group III-V, II- VI, and I- VII semiconductors us- 
ing the first-principles band structure method. We stud- 
ied the effect of the pi/ 2 local orbitals on the calculated 
Aso- The general trends of Aso of the semiconduc- 
tors are revealed and explained in terms of atomic SO 
splitting, volume deformation-induced charge renormal- 
ization, and cation-anion p-d couplings. In most cases, 



our calculated results are in good agreement with the ex- 
perimental data. The differences between our calculated 
value for a-Sn and HgTe, and to a lesser degree for InAs 
and GaSb, are highlighted. Experiments are called for to 
test our predictions. 
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